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Abstract 

The theory describing the far-field emission from a dipole source embedded inside a chiral sculptured 
thin film (CSTF), based on a spectral Green function formalism, was further developed to allow for 
infiltration of the void regions of the CSTF by a fluid. In doing so, the extended Bruggeman homoge- 
nization formalism — which accommodates constituent particles that are small compared to wavelength 
but not vanishingly small — was used to estimate the relative permittivity parameters of the infiltrated 
CSTF. For a numerical example, we found that left circularly polarized (LCP) light was preferentially 
emitted through one face of the CSTF while right circularly polarized (RCP) light was preferentially 
emitted through the opposite face, at wavelengths within the Bragg regime. The centre wavelength for 
the preferential emission of LCP/RCP light was red shifted as the refractive index of the infiltrating 
fluid increased from unity, and this red shift was accentuated when the size of the constituent particles 
in our homogenization model was increased. Also, the bandwidth of the preferential LCP/RCP emission 
regime decreased as the refractive index of the infiltrating fluid increased from unity. 
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1 Introduction 

Chiral sculptured thin films (CSTFs) constitute a remarkable class of engineered materials. A CSTF com- 
prises an array of parallel helical nanowires which can be grown on a substrate using vapour deposition 
techniques [1, 2]. By careful control of the fabrication process, both the optical properties and the multi- 
scale porosity of the CSTF can be tailored to order [3]. Accordingly, CSTFs are promising candidates as 
platforms for optical sensing, as well as a host of other applications [4, 5, 6]. 

Three quite different methods of utilizing CSTFs for optical sensing applications have been proposed. 
The first relies on the circular Bragg phenomenon, which CSTFs exhibit just as cholesteric liquid crystals 
do [7]. That is, within a wavelength regime known as the Bragg regime, a structurally right-handed CSTF 
almost completely reflects normally incident right circularly polarized (RCP) plane waves whereas normally- 
incident left circularly polarized (LCP) plane waves arc reflected very little. And vice versa for a structurally 
left-handed CSTF. The spectral shift in the circular Bragg phenomenon induced by infiltration of the void 
regions between the CSTF's nanowires may be exploited for sensing applications [8]. The second sensing 
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methodology is based on the excitation of surface-plasmon-polariton waves at the interface of an infiltrated 
CSTF and a suitable metal [9] . The third methodology for CSTF-based optical sensing — which provides the 
main motivation for the present communication — involves the emission of radiation from a chemiluminescent 
source embedded within a CSTF [10]. 

Chemiluminescent radiation may be generated within a CSTF — and harnessed for optical biosensing — 
as follows. Suppose that certain biorecognition sites are immobilized on the surface of a CSTF's nanowircs, 
and the CSTF becomes infiltrated with a solution containing complementary analytc biomolccules. The 
binding of the analytes to the biorecognition sites, perhaps in the presence of a transition metal com- 
plex, results in the generation of chemiluminescent photons [11, 12]. For example, the biorecognition sites 
and the analytes could be fragments of complementary single-stranded DNA, which combine to produce 
chemiluminescence in the presence of a ruthenium complex [13, 14, 15]. Inspired by the prospects of this 
chemiluminescent sensing scenario, the theory of emission from a dipole source embedded inside a CSTF 
was recently established [16]. In the following, we extend the theory in order to investigate the effects on 
infiltration upon the far-field dipole radiation. 

As regards notation, vectors and matrixes are represented in boldface, with the Cartesian unit vectors 
given by u x ,y,z\ dyadics are double-underlined; and both 4-vectors and 4x4 matrixes are enclosed within 
square brackets. The inverse of matrix [Z] is written as [Z] , while the inverse, adjoint and trace of the 
dyadic Z_ are written as Z_ _1 , Zl^ ancl tr (Z), respectively. The permittivity and the permeability of free 
space are denoted by e and [i , respectively; k = Lu^/e fi is the free-space wavenumbcr; A = 27r/fc 
is the free-space wavelength; and % = sj /i / e is the intrinsic impedance of free space. An exp(— iuit) 
time-dependence is implicit, with u being the angular frequency. 

2 Theory 

In this section we present an overview of the theory which describes the emission of electromagnetic radiation 
from a dipole source embedded inside a CSTF that is infiltrated by a fluid of refractive index m. The theory 
rests upon two recently-established pillars: First, in §2.1 we describe how the constitutive parameters of an 
infiltrated CSTF can be estimated using an inverse homogenization technique. Second, in §2.2 the far-field 
radiative emission from a dipole source embedded in an infiltrated CSTF is estimated via a spectral Green 
function approach. Comprehensive details of these two pillars arc available elsewhere [8, 16, 17]. Let us note 
that the theory presented in §2.1 represents an extension of that presented in [8] insofar as here the extended 
Bruggeman homogenization formalism is implemented. Also, the theory presented in §2.2 differs from that 
developed in [16] insofar as here we consider an infiltrated CSTF whose upper surface is immersed in a fluid 
of refractive index nt. A schematic diagram of the scenario under consideration is provided in Fig. 1. 

2.1 Constitutive parameters of an infiltrated CSTF 

A CSTF consists of an array of parallel nanowircs. The 'chirality' of the CSTF stems from the fact that the 
nanowircs are helical [1]. Such structures can be grown on a planar substrate — parallel to the plane z = 0, 
say — by the deposition of an evaporated bulk material. The helical shape of the nanowires is achieved by 
means of uniform rotation of the substrate about the z axis during the deposition process. We suppose that 
the deposited material is an isotropic dielectric material of refractive index n s . Notice that the value of n s 
needs to be determined as it may differ somewhat from the refractive index of the bulk material that was 
evaporated, depending upon the precise nature of the deposition conditions [18, 19, 20]. 

It is convenient here to regard the individual helixes of a CSTF as strings of highly elongated ellipsoidal 
inclusions, wound end-to-end around the z axis [21, 22]. The position vector 

r(d,Lp) = U_ ■ [ sin $ cos <p u„ + cos t9 u r + sin $ sin tp ] , $€[0,7r], ip € [0,2tt] (1) 

prescribes the surface of a particular ellipsoid relative to its centroid. Herein 77 is a linear measure of size 
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while the shape dyadic 

£ = u„ u„ + 7 T u T u r + 7 b u 6 u 6 

is expressed in terms of the normal, tangential, and binormal basis vectors per 

u„ = -u x sinx + u z cosx 
u r = u x cos x + u z sin \ 



(2) 



(3) 



with x denoting the inclination angle relative to the xy plane. An elongated ellipsoidal shape is achieved 
by selecting the shape parameters 76 > 1 and 7 T 3> 1. Since increasing j T beyond 10 does not result in 
significant effects for slim inclusions [22], the value 7 T = 15 is taken for the numerical results presented in 
§3. 

The proportion of a CSTF's total volume occupied by helical nanowires is represented by / £ (0,1). 
That is to say, the volume fraction 1 — / of a CSTF is not occupied by nanowires. 

Let us suppose now that the CSTF under consideration occupies the region — L < z < L, and is 
unbounded in extent in directions perpendicular to the z axis. At length scales much greater than the 
nanoscale, the CSTF is characterized by the relative permittivity dyadic 



£ +t = S 

— est J =z 



tt(z + L) 



n 



— y \/V./ — cg £j — y \/\,/ — , 



tt(z + L) 



z £ [-L, L] , 



with the rotation dyadics 



±i y (x) = u y Uy + (u x u x + u z u z ) cosx + (u z u x - u x u z ) sinx 
S_ (a) — u z u z + (u x u x + u y u y ) cos a + (u y u x — u x u y ) sincr 



(4) 



(5) 



The structural period is 2J7, and the handedness parameter h = +1(— 1) for a structurally right (left)-handed 
CSTF. The reference relative permittivity dyadic characterizes the local orthorhombic symmetry; i.e., 



-cstf 



' U n U n 



€bu U T U T + € cu Ub Ub , 



(6) 



with v = 1 denoting an uninfiltratcd CSTF (in which case the void regions between nanowires are assumed 
to be vacuous) and v — 2 denoting a CSTF in which the void regions are filled with a fluid of refractive 
index n£. 

We are required to estimate the relative permittivity parameters {e a 2, £62, £c2} for an infiltrated CSTF 
from a knowledge of the corresponding parameters {e a i, e&i, e c i} for an uninfiltratcd CSTF. A two-step 
strategy is employed. The first step is the estimation of the nanoscale parameters {n s , /, 7&} — which are 
not readily determined by experimental means — from a knowledge of {e Q i, £61, £ei}- As described in detail 
elsewhere [17], this can be achieved by applying the inverse Bruggeman homogenization formalism. Once 
{ n sj/i76} have been estimated, the second step can be taken wherein these parameters characterizing the 
uninfiltratcd CSTF are combined with {ti^,7 t } in order to determine the relative permittivity parameters 
{ea2, e&2j £02} for the infiltrated CSTF, by applying the Bruggeman homogenization formalism in its usual 
forward sense [8, 22]. For this second step we implement the extended version of the Bruggeman formal- 
ism which takes into account the nonzero size of the ellipsoidal inclusion particles [23]. That is, the size 
parameter r\ is taken to be small relative to wavelength(s) but nonzero. Details of the extended Bruggeman 
formalism, applicable to the locally anisotropic dielectric materials under consideration here, is provided in 
the Appendix. 



2.2 Radiation from dipole source inside an infiltrated CSTF 
2.2.1 Spectral Green function formulation 

Let us now introduce an electric dipole of moment p/2, oriented in the direction of the unit vector uj, which 
is embedded inside the CSTF at r = du z , d £ (-L, L). The corresponding source current density phasor is 
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given by 



J(r,u>) = — — UjS(z- d)S(x)6(y) , 



(7) 



where S(-) is the Dirac delta function. With a view to implementing a spectral-Green-function formalism 
[24] , the spatial Fourier transform representation 

r 2ir 



1 />OG rZTT 

J(r,oj) = — - / dun / dip j(z, ft, ip, ui) exp [in (x cosip + y sin ■0)1 
47r Jo Jo 



(8) 



is adopted, with j(z, n,ip,Lo) = — (iusp/2) uj 5 (z — d) in accordance with eq. (7). 
Our attention is focused on a single spatial-Fourier component of J(r,w), i.e., 

J(r, K, ip, ui) = j(z, K,ip,ui) exp [in (x cosip + y simp)] , z 6 (— L, L 

Similarly, we write the electromagnetic field phasors inside the CSTF as 

E(r, ft, ip, uj) ~ e(z, K,ip,uj) exp [zk (x cos ip + y sin -0)] , "I 

H(r, n,ip,uj) = h(z, k, ip, ui) exp [in (x cos ip + y sin ?/>)] J 

The procedure whereby the Fourier components (9) and (10) are combined with the constitutive relations 
for the CSTF and the frequency-domain Maxwell curl postulates, and then solved to find expressions for 
the phasors inside the CSTF, is comprehensively described elsewhere [16]. Therefore, here we simply state 
the particular solution 



z e (-L, L) 



(9) 



(10) 



[f(L,K,V,cj)] = [M(L ) K,V,w)][f(-L,K,V,w)] 



[M(L, k, ip, uj)] [M(d, k, ip, uj) 



^ [g{d,K,ip,u;)] 



wherein the column vectors 



[{(z,K,ip,u)] = 



e(z,K,tp,u;) ■ u x 
e(z, K,ip,w) ■ u y 
h(z, K,tp,u) ■ u x 

h(Z,K,tp,Uj) ■ Uy 



(11) 



(12) 



and 



i(z,K,ip,u>)] = 



cos 2 x + tbv sin x 



K COS Ip 

k sin ip 



(e i/ - tbv) sinx cosx sin 
(e^y - t av ) sinx cosx cos 



Mz+L) 

' n 
1 n 






Uj Uy 

-uj • u x 



(13) 



The expressions for the 4x4 matrizants \M.(L, n,ip,ui)] and [M(d, ft, ip, uj)] — which are straightforwardly 
derived, but too cumbersome to reproduce here — are available in standard works [1, Chap. 9]. We note 
that the piecewise uniform approximation technique provides a convenient method for their evaluation [1, 
Chap. 9]. 



2.2.2 Boundary value problem 

Next we turn to the two half-spaces z < —L and z > L. The half-space z < — L is vacuous while the 
half-space z > L is filled by a fluid of refractive index rig. In consonance with cqs. (10), the electromagnetic 
field phasors for these two half-spaces may be expressed as 

E(r,K,tp,ui) = —= \-b L (K,ip,ui) (is - p_) + b R (n,ip,uj) (is + p_)] 



U(r,K,ip,u) 



V2 

x exp{i [k (x cosip + y simp) — a (z + L)]} 

[b L (ft, ip, uj) (is - p_) + b R (ft, ip, ui) (is + p_)] 



x exp{z [k (x cos ip + y simp) — a (z + L)]} 



z < -L, 



(14) 
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and 



E(r,/e,V,w) = -j= [c L (K,ip,w) (is - p+) - c R (n,ip,uj) (is + p+)] 



x exp {i [k (x cos ip + y sin?/;) + (z — L)}} 



U(T,K,1p,u) = — 



[c L (k, ip, uj) (is - p+) + c R (k, ip, uj) (is + p+)] 



mV2 

x exp {i [k (x cos ip + y sin?/;) + at (z — L)}} 



z > L, 



(15) 



where a = +y/k% — k' 2 , = +yhj — n 2 , fc^ = fc n^ and 77^ = r] /ri(. The complex-valued amplitudes 
fe^ (re, ip, uj) and (k, ip, uj) represent the LCP components, while bn (k, ip, uj) and cr (k, ip, uj) likewise rep- 
resent the RCP components. The unit vectors 



and 



P± = 



-Us; sin ip + u y cos ip 



T(at /k ) (u x cosip + u y sin?/;) + (n/k ) u z , 



< -L 



T{a e /h) (u x cos ip + u y sin ip) + (n/ke)u z , z > L 



(16) 
(17) 



relate to the perpendicular- and parallel-polarization states of the plane wave, respectively. 

The tangential components of E(r, k, ip,oS) and H(r, K,ip, uj) are required to be continuous across the 
pupils at z = — L and z = L of the CSTF. Thus we have that 



and 



[f(-L,K,ip,uj)] = -j= [K(K,ip,uj,a ,k ,r) )] 



[i(L,K,ip,uj)} = — [K(K,ijj,uj,a e ,ki 
v2 






—i [b L (k, ip, uj) - b R (k, ip, uj)} 
b L (n,ip,u>) + b R (n,ip,uj) 

i [c L («, ip, uj) ~ cr («, ip, uj)} 
- [cl (k, ip, uj) + c R (k, ip, uj)} 







(18) 



(19) 



where 



— sin^ 
cos^ 

_ Q COS -0 

kr} 
_ a sin ip 
kr} 



a cos ip 

a sin ip 
k 

sin y: 
cos ip 



— sini/j 
cosifj 



a cos -ip 
k 

a sin w 
k 

a cos ip sin ip 

k-7] Tj 

a sin ip cos ip 

k7] 77 



(20) 



The unknown four amplitudes b^^R (n,ip,uj) and cl.r (n,ip,uj) can now be determined using standard alge- 
braic manipulations, by combining eqs. (18) and (19) with the particular solution (11). 



2.2.3 Emitted far field phasors 

The emitted electromagnetic phasors in the half-spaces z < —L and z > L are found by summing the 
corresponding spatial Fourier components (14) and (15) per 



E(r,w) = — 77 / dn / dip K E(r, n, ip, uj) 
4tH Jo Jo 

H(r,w) = — - / dn I dip «H(r, K,ip,uj) 

4tt 2 Jo Jo > 



z i [-L, L]. 



(21) 
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Asymptotic approximations to these integrals, representing the emitted field phasors in the far zone, are 
provided by [25, 26] 

icos6 obs r .„»../. „>,. „^„n .„>,./. „>,. j,.m exp (ik r°} s ) 



and 



E(r obs ,w) ft 
H(r obs ,w) r 

E(r obs ,u;) 
H(r obs ,u;) 



2y/2n 
cos8 obs 

icos9 obs 
2^2-k 



[-bf s (is obs - p obs ) + b°£ s (is obs + p° bs )] 



fe f° f 



~obs 



k f°} s 



[bf s (is obs - p° bs ) + bf s (is obs + p° 6s )] exp ( zfe ° r - ) 
[cf 8 (is obs - p obs ) - (is obs + p ob s )] 



z ot,s < _ L ^2) 



exp (iktr°l 



^obs f -„obs n obs\ I „obs (•„obs , n obs 



pf)] 



k t f obs 
exp (ifarf 8 ) 



z obs > L, (23) 



wherein r± = |r + Lu z | and the superscript obs denotes evaluation at the distant observation point 



obs — ( 'obs nobs n l,ob 



e obs ,ip obs ) with z obs = r obs ■ u 2 . Note that the approximations (22) and (23) are appropriate at 
distances far from the CSTF pupils but not in the vicinity of obs = ir/2 [25, 26]. 

For practical purposes, the radiation field in the far zone is conveniently characterized in terms of the 
time-averaged Poynting vector, which we write as the sum of LCP and RCP contributions per 

P(r obs , u) = P iC p(r ob '» + P RC p(r obs , w), 



wherein 



P LC p(r obs ,^) ft { 



1 

2^o 



( cos e ob 



\2irk " ob 



j.obs 



1 



( cos 9 



obs \ 2 



and 



P^Cp(r obs , W )ft i 



1 



ir(K s \ 2 ) 



i 



obs 1 2 s ! 



I 2 H 



\2irk e rf s 



cos a 
2nk r°_ bs 

cos6 obs 
2nkff obs 



'lobs 



z obs < _ L 



z obs > L 



(24) 



(25) 



^obs „obs 



'xobs 



z obs < _ L 



z obs > L 



(26) 



at the observation point r° 



3 Numerical results 

For our numerical investigations we selected a structurally right-handed CSTF (i.e., h = +1) with struc- 
tural half-period f2 = 200 nm. Experimentally-determined values for the relative permittivity parameters 
{fai, e&i, £ci}) which characterize the uninfiltrated scenario, were used. In light of the absence of appropriate 
data for CSTFs, we chose the relative permittivity parameters 



Col 



£61 



e c i 



1.0443 + 2.7394 



1.6765 + 1.5649 



1.3586 + 2.1109 



2x^ 

7T 
7T 
IT 



1.3697 



0.7825 



1.0554 



2Xv 

7T 



(27) 
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with 

tan x = 2.8818 tan Xv , (28) 

which were determined by measurements on a columnar thin film made from patinal® titanium oxide 
[27, 28]. The angle Xv (radians) in eqs. (27) represents the average direction of the vapour flux relative to 
the substrate during the deposition process. As described in a previous study [17], the inverse Bruggeman 
homogenization formalism yields the corresponding nanoscale model parameter values: n s = 3.0517, / = 
0.5039 and jb = 1-8381 for the vapour flux angle Xv = 30°. 

The relative permittivity parameters {e Q 2, £b2, £c2} for the infiltrated CSTF, as computed using the ex- 
tended Bruggeman homogenization formalism, arc graphed as functions of the refractive index of the in- 
filtrating fluid rn G (1.0,1.5) and the relative size parameter k a r\ £ (0,0.2) in Fig. 2. The real parts of 
{£q2, £62, £02} increase uniformly as ng increases; they also increase uniformly as rj increases but more slowly. 
The imaginary parts of {e a 2 7 £b2, £02} are fairly insensitive to ng, but these quantities increase exponentially as 
the size parameter increases from zero. The manifestation of constitutive parameters with nonzero imaginary 
parts for homogenized composite materials, when the component materials are themselves nondissipative, is 
a well-recognized phenomenon. This phenomenon — which arises in higher-order homogenization theories, 
such as the strong-property-fluctuation theory [29, 30] and extended variants of the Bruggeman and Maxwell 
Garnett formalisms [31, 32, 33] — may be attributed to radiative scattering loss from the macroscopic co- 
herent field [34]. In particular, we note that Ime„2.62.c2 — > in the limit as 77 — > 0. 

Now let us turn to the radiation emitted from a dipole source embedded within the CSTF. We chose 
the source to be located relatively close to the upper surface of the CSTF at d = L — 40 nm; and the 
source orientation was given by uj = S_ (d, h) ■ u„ 2 . From the point of view of sensor applications, we are 
particularly interested in the ability of the CSTF to discriminate between LCP and RCP light. Accordingly, 
the wavelength ranges considered here were selected to include the circular Bragg regime. The location and 
extent of the the circular Bragg regime were conveniently estimated by the centre wavelength [1] 



Af {e av ,e hv ,e cv ,6 obs ) ~ Q \^A^\ + 
and the full-width-at-half-maximum bandwidth [1] 



^av^-bv 



cos 2 x + tav sin 2 x 



\ / cos6 obs , 



(29) 



(AAof {e av ,e bv , e c „, 9 obs ) ~ 2Q v^~| - 



cos 2 x + ^av sin 2 x 



(30) 



The thickness ratio L/il = 30 was chosen to ensure that the circular Bragg phenomenon is fully developed. 

In order to appreciate the effects of infiltration, we must first consider the uninfiltrated scenario. In Fig. 3, 
the projections of [P Lcp{r° bs , oj) I and \P Rcp( r ° bs , w )| onto the z = plane are mapped for rig = 1 (i.e., 



1). The contributions to the time-averaged Poynting vector were scaled by a factor of 10 



13, 



IpI" 



The radiation emitted through the CSTF's pupils at z — L and z = —L are both represented. Results are 
presented for A = Xf r (e a i, e&i, e c i, 0°) = 763.6 nm, as well as for A = Xf r (e a i, e^i, e c i, 0°) — 80 nm and 
A = A,f r (e a i, efei, e c i, 0°) +40 nm. Also presented in Fig. 3 are evaluations of the real-valued parameter 



10 



16 



d6° 



dij° 



„obs 



nobs 



^obs = Q 



|P,(r obs , W )| , j e {LCP, RCP} 



(31) 



wherein p x = 0, p 2 = 0.95tt/2 for z obs > L, and p 1 =ir- (0.95tt/2), p 2 = n for z obs < -L. The quantity Tj 
delivers a measure of the total rate of energy flow into the half-spaces z > L and z < —L. 

There are several notable features in Fig. 3, especially concerning differences between the LCP and 
RCP emission characteristics, which relate to the circular Bragg phenomenon: (i) For emission which is 



2 We also considered the source orientations uj = S (d, h) ■ u r and uj = S (d, h) ■ 115. The results for these cases are not 
presented here as they were qualitatively similar to the case uj = S_ (d, h) ■ u n . 



7 



approximately normal to the two pupils of the CSTF, RCP radiation is preferentially emitted through the 
pupil at z — L for A = 763.6 nm whereas at the same wavelength LCP radiation is preferentially emitted 
through the pupil at z = —L. (ii) At A = 683.6 nm, the differences between LCP and RCP emission are 
very small for | cos9 obs \ ~ 1, but we can see that RCP radiation is preferentially emitted through the upper 
pupil whereas LCP radiation is preferentially emitted through the lower pupil for 0.4 < | cos8 obs \ < 0.8. This 
observation is in accordance with the blue shift of the circular Bragg phenomenon for obliquely incident plane 
waves [1]. (iii) The distinction between the LCP and RCP patterns is barely noticeable when | cos o b s \ < 0.4, 
for all three wavelengths considered. This is indicative of the severe diminishmcnt of the circular Bragg 
phenomenon for highly oblique planewave incidence [35]. (iv) There is very little evidence of the CSTF 
discriminating between LCP and RCP radiation at A = 803.6 nm; this wavelength represents an upper 
bound on the circular Bragg phenomenon for all angles of incidence. 

Let us remark too upon the distinctive pattern of concentric rings that appears in Fig. 3. These are 
Fabry-Perot interference rings, arising due to the finite thickness of the CSTF. Indeed, even if the CSTF 
were replaced by a homogeneous isotropic dielectric material the ring pattern would still be observed [16]. 

The results of infiltration can be observed in Figs. 4 and 5, which show the projections of \PLCp{ r ° bs , w ) | 
and |Pflcp( r0bS j w )| (as before, scaled by a factor of 10 13 w -2 ) onto the z — plane for ng = 1.25 
and ri£ — 1.5, respectively. The relative permittivity parameters {e a 2, e&2, £c2} for the infiltrated CSTF were 
computed using the non-extended version of the Bruggeman homogenization formalism (or, equivalcntly, the 
extended version with rj = 0). In keeping with Fig. 3, results are presented for A = A,f r (e a 2, £b2, £c2, 0°), as 
well as for A = A<f r (e a 2, £b2,£c2, 0°) — 80 nm and A = A,f r (e a 2, £62, £c2, 0°) + 40 nm. As for the uninfiltrated 
case, we see that RCP radiation is preferentially emitted through the pupil at z = L at whereas LCP radiation 
is preferentially emitted through the pupil at z = —L, at the centre Bragg wavelength A<f r (e Q 2, £b2, £ C 2, 9 obs ) 
for | cos8 obs \ ~ 1. Two effects in particular of infiltration are apparent from Figs. 4 and 5 as the refractive 
index of the infiltrating fluid increases: (i) the distinction between LCP/RCP emission for |cos0 obs | ~ 1 
becomes increasingly red shifted; and (ii) the blue shift in the distinction between LCP/RCP emission for 
oblique angles of incidence becomes less pronounced. 

In fact, the main effects of infiltration are quite well predicted via the empirical relations (29) and (30). 
To see this, the centre Bragg wavelength and the full-width-at-half-maximum bandwidth arc plotted versus 
ri£ in Fig. 6. As in Figs. 4 and 5, the non-extended version of the Bruggeman homogenization formalism was 
used to estimate the relative permittivity parameters {e a 2, £(,2, £02}- The red shift in the centre wavelength, 
and the reduction in the full-width-at-half-maximum bandwidth, resulting from the refractive index of the 
infiltrating fluid being increased are obvious from Fig. 6 — and these effects are in complete agreement with 
the emission patterns observed in Figs. 3-5. 

We consider now the influence of the linear size of the ellipsoidal particles which represent the CSTF's 
helical nanowires, per the homogenization model described in §2.1. As a representative example, we repeat 
the calculations of Fig. 4 but here using the extended Bruggeman homogenization formalism to estimate 
the relative permittivity parameters {e a 2, Cb2, £02} of the infiltrated CSTF. We set the size parameter to 
be 77 = 0.1/fc o . The corresponding projections of |PLCp(r obs ,u;)| and \Pjfcp{ r ° bs : w )| (as before, scaled 
by a factor of 10 13 w -2 \p\ 2 ) onto the z = plane are provided in Fig. 7 for rig = 1.25. By comparing 
Figs. 4 and 7, we see that the size parameter has a relatively minor but significant influence on the emission 
patterns. The centre wavelength for the distinction between LCP/RCP emission is slightly higher in Fig. 7 
as compared to Fig. 4. Furthermore, the total energy flux emitted from the CSTF — as estimated by the 
scalar parameter r lcp,rcp — is substantially smaller when we consider 77 = 0.1/A: as opposed to r\ = 0. 
This is most noticeable for radiation emitted through the CSTF pupil at z = —L, which is a consequence of 
the dipole source being closer to the z = L pupil. 

The main effects of the size parameter 77 may be estimated quite well using the empirical relations (29) and 
(30). In Fig. 8, Af r (e a2 , e&2, e c2 , 0°) and AA Br (e a2 , e h2 , e c2 , 0°) arc plotted versus k rj for nt S {f -0, 1.25, 1.5}. 
The modest increase in the centre Bragg wavelength is clear; and we note that the rate of increase is greatest 
when the refractive index of the infiltrating fluid is smallest. The full-width-at-half-maximum bandwidth 
is relatively insensitive to the size parameter 77, regardless of the value of ri£. Qualitatively similar results 
are found when angles 8 obs > 0° are considered. 
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4 Closing remarks 



Using a spectral Green function formalism, in conjunction with an inverse homogcnization formalism, the 
effect of infiltration on the emission from a dipole source embedded within a CSTF has been characterized. 
Based on numerical studies, our conclusions may be summarized as: 

• The centre wavelength for the preferential emission of LCP/RCP radiation is red shifted as the refrac- 
tive index of the infiltrating fluid increases from unity. Furthermore, the red shift is accentuated when 
the size of the ellipsoidal particles which represent the helical nanowircs of the CSTF is increased. 

• The bandwidth of the preferential LCP/RCP emission regime decreases as the refractive index of the 
infiltrating fluid increases from unity. 

• The main effects of infiltration may be reasonably predicted using the simple empirical formulas (29) 
and (30) which provide estimates of the centre Bragg wavelength and the corresponding full-width- 
at-half-maximum bandwidth. 

Through the elucidation of the effects of infiltration, a further step towards the practical realization of 
CSTF-based optical sensors, including biosensors which harness chcmiluminesccnt emission from inside a 
CSTF, has been taken. 
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Appendix 

In order to estimate the reference relative permittivity dyadic of an infiltrated CSTF, a local homogcnization 
procedure is carried out within a plane parallel to z = 0. Two components are to be homogenized: (i) a 
planar array of similarly-aligned ellipsoidal particles — of refractive index n s and specified by the linear 
size parameter rj and shape dyadic U_ per eq. (1) — which are the building blocks of the helical nanowircs; 
and (ii) void regions — represented as a collection of spherical particles, of radius specified the linear size 
parameter rj — which are infiltrated with a fluid of refractive index ng . The linear size parameter 77 is taken 
to be much smaller than the wavelength(s) but not vanishingly small. 
From the nonlinear Bruggeman equation [23] 



(i-/)g, /(2) +/g s/(2) = a, 

(2) 

e cst y can be extracted by standard numerical methods, such as the Jacobi technique 
polarizability density dyadics 



(32) 

Herein, the 



a/zoo = W^l-iclf 



i s /( 2 ) = [V^L~i% 



(2) 
cstf 



> - 



(33) 



where the depolarization dyadics 



1 
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dip I d$ sin 0(W° + ifW + ^ 



(34) 
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with 



A = sin 2 -& cos 2 ip u„ u„ 



cos 2 "D 



sin 2 $ sin 2 ip 



> . 



(35) 



-A = sin # cos 2 tp u„ u„ + cos 2 ■& u T u T + sin -dsm (pUbUb 
The dyadic integrands in eq. (34) are given as 
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and scalar quantities 
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(40) 
(41) 
(34) [37]. 



Numerical methods are generally needed to evaluate the surface integral on the right side of eq. 
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Figure 1: A schematic diagram of the scenario under investigation: radiation emitted from a dipole source 
embedded within an infiltrated CSTF occupying —L < z < L. 



13 




(e) (f) 

Figure 2: Real and imaginary parts of the relative permittivity parameters e a 2, tb2 and e C 2 versus refractive 
index ng G (1, 1.5) and relative size parameter k rj £ (0, 0.2). 
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Ao = 683.6 nm r LCP = 26.7 Ao = 763.6 nm r LCP = 32.3 Ao = 803.6 nm r LCP = 37.6 




(j) (k) (1) 

Figure 3: Projections of |Plcp| an d |P.rcp| (scaled by 10 13 u;~ 2 |p|~ 2 ) onto the z = plane for z obs > 
and z obs < —L. Here A G {683. 6nm, 763. 6nm, 803. 6nm} and ni = 1.0. 
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A = 745.2 nm r LCP = 14.9 A = 825.2 nm r LCP = 18.9 A = 865.2 mil l',, n , = 21.3 




(j) (k) (1) 

Figure 4: As Fig. 3 but with A € {745.2nm, 825. 2nm, 865.2nm} and n; = 1.25. The relative permittivity 
parameters {e a 2, e&2, £02} for the infiltrated CSTF were computed using the non-extended version of the 
Bruggeman homogenization formalism. 
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A = 804.9 nm r LCP = 10.7 A = 884.9 mil r LCP = 13.3 A = 924.9 nm r LCP = 15.0 




(j) 00 (i) 

Figure 5: As Fig. 4 but with A e {804.9nm, 884.9nm, 924.9nm} and n ; = 1.5. 
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Figure 6: X^ r and AA^ r plotted against rc z for 9 obs = 0° (solid, red curve), 6 obs = 30° (dashed, green 
curve and e obs = 60° (broken dashed, blue curve). The relative permittivity parameters {e a 2, £b2, £c2} for 
the infiltrated CSTF were computed using the non-extended version of the Bruggeman homogenization 
formalism. 
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A = 746.5 nm r LCP = 13.9 



A = 826.5 nm r LCP = 17.1 



A = 866.5 nm r LCP = 19.2 
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Ao = 746.5 nm r RCP = 14.7 
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Figure 7: As Fig. 4 but the relative permittivity parameters {e a 2, £62) e c2} for the inhltrated CSTF were 
computed using the extended version of the Bruggeman homogenization formalism with rj = 0.1/fc o , and 
A € {746.5nm,826.5nm,866.5nm}. 
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Figure 8: Xq t and AXg r plotted against k rj for ri( = 1.0 (solid, red curve), ri( = 1.25 (dashed, green curve) 
and n,( = 1.5 (broken dashed, blue curve), as computed using the extended Bruggeman homogenization 
formalism. The angle 6 obs = 0°. 
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